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Abstract 

This paper is devoted to integrability conditions for systems of linear difference and differential 
equations with difference parameters. It is shown that such a system is difference isomonodromic 
if and only if it is difference isomonodromic with respect to each parameter separately. Due to 
this result, it is no longer necessary to solve non-linear difference equations to verify isomon- 
odromicity, which will improve efficiency of computation with these systems. 
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1. Introduction 

In this paper, we improve the algorithm that verifies if a system of linear difference and 
differential equations with difference parameters is isomonodromic (Definition 3.1). Given a 
system of difference or differential equations with parameters, it is natural to ask if its solutions 
satisfy extra linear difference equations with respect to the parameters. Consider an algorithm 
whose input consists of a field K with commuting automorphisms <j>\, . . . , <f> q , cr\, . . . , o~ r and 
derivations d\,...,d m onK and a system of difference equations 

< /> 1 (Y)=A 1 Y,..., ( /, q (Y)=A q Y, d\Y = B\Y,..., d m Y = B m Y, (1) 

where A, e Gh n (K) and B ; e M n (K), 1 ^ i ^ q, 1 ^ j ^ m. Its output consists of such 
additional linear difference equations 

o-iY = CiY,...,o- r Y = C r Y, (2) 

where C, e Gh„(K), 1 ^ i ^ r, if they exist, for which there exists an invertible matrix solution 
of (1) that satisfies (2) as well. For an invertible matrix solution of (1) to possibly exist, the A,'s 
and B /& must satisfy the integrability conditions (5), (6), and (7). Moreover, existence of the 
C;'s above is equivalent to the existence of C,'s satisfying another large collection of integrabil- 
ity conditions (9), (10) and (11). Such systems (1) are called isomonodromic in analogy with 
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differential equations [31, 32, 23, 24, 19]. Isomonodromy problems for g-difference equations 
and their relations with g-difference Painleve equations were studied in [25, 26]. 

The main result, Theorem 4.1, states that (11), which are non-linear difference equations 
in the C,'s, do not have to be verified to check the existence of the C,'s. More precisely, the 
existence of C,'s satisfying (9) and (10) implies the existence of new invertible matrices that 
satisfy all of (9), (10), and (11). 

Since, due to our result, we only need to check the existence of solutions (that have entries 
in the ground field K) of a system of linear difference and differential equations, a complexity 
estimate for verifying whether a system of difference and differential equations is isomonodromic 
becomes possible due to [1, 2, 4, 6, 7, 8]. 

A similar problem but for systems of differential equations was considered in [ 19], motivated 
by the classical results [23, 24]. Differential categories developed in [18] formed the main tech- 
nical tool in [19]. On the contrary, our proofs are constructive, which makes them more suitable 
for practical use. In particular, from our proof, one can derive an algorithm that, given a com- 
mon invertible solution of (9) and (10), computes a common invertible solution of all (9), (10), 
and (11). 

If one attempts to find a common solution for all of (9), (10), and (11) in a naive way, one 
might have to adjoin new elements to the base field K that are not constant with respect to the 
derivations and automorphisms, which is not desirable. On the other hand, Theorem 4. 1 shows 
how one can limit, in a computable way, these newly adjoint elements to those that are constant 
with respect to all derivations and the automorphisms that do not correspond to the parameters, 
which is another advantage of this approach. This is also different from the approach taken in 
[5] to give an explicit treatment of a similar problem but for differential parameters, where a 
linearly differentially closed assumption is imposed on the field of constants with respect to the 
automorphisms. 

The main result is expected to have further applications. In parameterized differential Galois 
theory [9], there are several algorithms for computing Galois groups. For 2x2 systems, they are 
given in [3, 17]. An algorithm, more general in terms of the order of the system, [28, 27] uses 
the differential analogues [19] of our results to make the computation more efficient. Our results 
may prove useful in the design of an algorithm for computing parameterized difference Galois 
groups with difference parameters, as it has happened in the differential case, once [16] or [30] 
is generalized to the case of several difference parameters, which is a challenge on its own. 

These and other difference and differential Galois theories [22, 21, 12, 14, 13, 15, 11, 10, 20, 
29, 33] have been developed to compute all difference and differential algebraic relations that 
solutions of systems of linear difference and differential equations satisfy and have found many 
practical applications. 

The paper is organized as follows. Notation is introduced and the basic notions of differential 
and difference algebra are reviewed in §2. The integrability conditions are discussed in §3, where 
the main result is illustrated by showing a concrete example, in which we perform computation 
with the basic hypergeometric ^-difference equation with parameters (Example 3.2). The main 
result and its proof are described in §4. 
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2. Basic definitions 



A A-ring is a commutative associative ring with unit 1 together with a set A = {<9i, . . . , d m } of 
commuting derivations 3,- : R — » such that 

3,(a + fo) = 3,(a) + di(b), di(ab) = di(a)b + adi(b), a,b e R, 1 < i < m. 

For example, Q is a {3}-field with the unique zero derivation. For every / € C(x), there exists a 
unique derivation d : C(x) — > C(x) with <9(x) = /, turning C(x) into a {<9}-field. 

For a ring R, M n (R) denotes the set of n x n matrices with entries in R, and GL n (R) are the 
invertible matrices in M„(R). 

A <l)-ring R is a commutative associative ring with unit 1 and a set O = {<j>\, . . . ,(f> q ] of 
commuting automorphisms tf>, ■ : R — » R, 1 ^ i ^ q. A \<t>, A}-ring is a <5-ring and a A-ring such 
that, for all e O and d e A, <pd = d<f>. For subsets <5' c <5 and A' c A, let 

R*'- A ' -={reR\ 0(r) = r, 5{r) = 0, for all e <&', 6 e A'}, 

the set of {O', A'}-constants of 

To discuss difference parameters in systems of difference and differential equations, it will 
be convenient for us to split the set of automorphisms into two subsets. In such a case, we will 
write {<!> U 2, A) to emphasize that the set <t> U 2 of given automorphisms is split into two subsets, 
O and 2. In this case, the 2 corresponds to the parameters. 

Let 2 = [ct\ , . . . , o>}, R be a 2-ring, and yi, . . . , y n be indeterminates over The ring 

R{y u . . . ,y„} E := * [fy |0 = < ■...■ofr.KK n, ij eZ,l^j^r] 
is called the ring of ^-polynomials, which is naturally a 2-ring. 

Example 2.1. To clarify the above, we will show a few basic examples of {<t> U 2, A}-rings: 

1. R = Q(xi,...,x m ,a 1 ,...,a r )with<l) = {<p\, . . . , (f> m ], 2 = {cr,, 07}, and A = {di,...,d m } 
defined by 

4>iif){x\, . . . ,x m ,a\, . . . ,a r ) = f(xi,...,Xi+ 1, . . . ,x m ,a u . . . ,a r ), f e R, U i < ra, 
o-i(f)(x u ...,x m ,a u ...,a r ) = f(x u ...,x m ,a 1 ,...,a i + l,...,a r ), f e R, U i ^ r, 
3; = <9/<9x;, 1 < i ^ m, 

where xi , . . . , x„„ ai, . . . , a r are transcendental over Q. 

2. # = Q(x, a) with <D = {0}, 2 = {cr}, and A = \d) defined by 

(f>(f)(x, a) = f(qix, a), cr(f)(x, a) = f(x, q 2 a), d = xd/d x , 
where + q\, q 2 e Q and x and a are transcendental over Q. 

Example 2.2. Difference and differential equations with difference parameters that appear in 
practice include, among many others, 

y"(x) + (x 3 + ax + b) ■ y(x) = 

and 

2 (a + b)x-(l +c/q) x-1 

abx - c/q abz - c/q 
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3. Integrability conditions 

In this section, the integrability conditions are introduced and an example illustrating our ap- 
proach is given. The main result will be shown in §4. Let K be a {<£ U X, A}-field of characteristic 
zero, where <B = {(pi, . . . , <f> q }, S = {cr\ , . . . , cr r }, and A = {<9i , . . . , d m ), and 

A\,...,A q € Gh n (K) and B u . . . ,B m e M„(2T). 

Consider the system of difference-differential equations 

1 (F)=A 1 7,...,0 ? (F)=A^ 5i(F)=Bi,...,a m (F) = J B m K (4) 

If L D K is a {O, A}-ring extension and Z € GL„(L) satisfies (4) then, for all i, j, 1 ^ i, j ^ g, we 
have 

4>j(AdAjZ = ^(0,(Z)) = = 0,(A ; )A,Z. 

Therefore, we obtain 

<t>j(Ai)Aj = 0,(A ; )A,- (5) 

Moreover, 

<f>j(BdAjZ = <Pj(di(Z)) = d0j{Z)) = di(Aj)Z + AjBiZ 

and 

<9,(fi,)Z + fi^Z = <9/di(Z)) = 5;(3/Z)) = di(Bj)Z + BjBjZ, 

which imply that 

(pj(Bj)Aj = di(Aj) + AjB h 1 < i < m, 1 < j < q, (6) 

and 

a/Bj) - di(Bj) = [B h Bi] , 1 < i, j < m. (7) 
If, in addition, there exist Ci, . . . ,C r e GL„(/T) such that 

CTi (Z) = CiZ, . . . , cr r (Z) = C r Z, (8) 

then, similarly to the above, 

<pj(Ci)Aj = ^(A^Q, 1 s$ i sC r, 1 j s$ q, (9) 

o-j(Bi)Cj = di(Cj) + CjBj, 1 < i < m, 1 < j < r, (10) 

tr/CflC, = criCC^C;, 1 sC /, j < r, (11) 

Therefore, (5), (6), (7), (9), (10), and (11) are necessary conditions for the existence of a common 
invertible matrix solution of (4) and (8) with entries in a field L. They are also sufficient. Indeed, 
let 

L = K(xu, ■ ■ - ,x nn ), 

with 

di((x rs )) = Bj(x rs ), 1 s$ i m, ; ((x„)) = A ; (x„), 1 < j < q, 
cr k {{x rs )) = C k (x rs ), 1 < k ^ r. 

Then (5), (6), (7), (9), (10), and (1 1) imply that L is a {€> U I, A}-field. 
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Definition 3.1. The system of linear difference-differential equations 

= A X Y, . . ., 4> q {Y) - A q Y, 8 1 (Y) = B 1 ,..., d m (Y) = B m Y 



with A x , . . . ,A q e GL„(K) and B u 
isomonodromic if there exist C\ , . . 



.,B m e M„(K) satisfying (5), (6) and (7) is called difference 
, C r e GL„(K) satisfying (9), (10), and (11). 



Connections to analytic interpretations of a similar notion, differential isomonodromy, and 
parameterized Galois theory can be found, for example, in [9, §5] and [19, §6.2]. 

Example 3.2. Consider the field K :- Q(x,a,b,c) as a O U S-field, where €> 
{cr a , o-f,, cr c ) and, for all / e K, 



{<Pa 



<p q (f){x, a, b, c) - f{qx, a, b, c), 
o~b(f)(x, a, b, c) - f(x, a, qb, c), 

For the equation (3), the matrices 



o~ a (f)(x, a, b, c) = f(x, qa, b, c), 
cr c (f)(x, a, b, c) = f(x, a, b, qc). 



Ci := 



1 -a 
aq(x — 1) aq(l — ax) + c(a — 1) 



abqx ■ 



abqx ■ 



C 2 :-- 



and 



C 3 := 



x(c(a + b) - ab) - c 

cx 
x-l 



1 -b 

bq(x - 1) bq{\ -bx) + c(b- 1) 
abqx — c abqx — c 

c — abx\ 



x 

c — abx 



X CX 

which can be calculated using the Rational Solution tool in the QDif f erenceEquations 
package of Maple (see also [30, Ex. 2.35] for the case a — b and c = q), satisfy (9). How- 
ever, (11) is satisfied only for the pair C\ and So, the ^-difference equation (3) is difference 
isomonodromic over any <t> U {cr a , <Ti}-field containing K. Applied to this situation, the proof of 
Theorem 4. 1 contains an algorithm that turns (3) into a difference isomonodromic equation with 
respect to <1> U 2 just by extending K by at most one cr c -transcendental element x\ to K, which is 
constant with respect to <\> q . In particular, a calculation in Maple using the above package shows 
that the solution space of (19) has dimension d — 1 and 

b — c 



O-a(xi) 



Since 



qa 



-xi, cr b (xi) 



1 



qb 



-Xi. 



(12) 



(a - c)(b — c) 

satisfies (12), following the proof of Theorem 4.1, we replace C3 by 

' x(c(a + b) — ab) — c 1 c — abx 



D* := 



(a 



■ c)(b — c)cx 
x-l 



(a - c)(b — c)x 
c — abx 



(a - c)(b — c)x (a - c)(b — c)cx) 

Thus, equation (3) is difference isomonodromic with respect to X over K with the matrices Cj, 
C2, and D3, and the proof of Theorem 4. 1 has helped us discover the latter matrix. 



Let L\ be a {O U S, A}-field and L2 be a S-field containing L, ' , then R 



,«,A L2 IS 



an integral domain [22, Lem. 6.11] and, therefore, L1L2 := Quot(/?) is a {O U 2, A}-field with 
{<!), A}-constants equal L 2 . 
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4. Main result 



This section contains the main result, Theorem 4.1, which allows one to reduce the number 
of integrability conditions to be tested. 

Theorem 4.1. Let A u . . . , A q e GL n (K) and B u ...,B m e M„(K) satisfy (5), (6), and (7). If 
there exist C\,...,C T e GL„(X) satisfying (9) and (10), then there exist 

1) a computable Y.-field F generated over K® ,& by at most (r — l)n 2 elements and 

2) Di, . . . ,D r e GL,„(KF) such that all integrability conditions are satisfied: 

<t>j(Di)Aj = o-i{Aj)Du 1 < i < r, l^j^q, (13) 

o-j(Bi)Dj = di(Dj) + DjBi, 1 < i < m, 1 ^ j < r, (14) 

°~ j(Di)Dj = o-i{Dj)Di, 1 < i, j sC r. (15) 

Proof. This will be done by induction. Let there exist C\,...,C r e Gh n (K) satisfying (9) 
and (10) and let k be such that 2 ^ k ^ r. Suppose we have already computed a S-field F k -\ 
generated over Zf°' A by at most (k - 2)n 2 elements and 

D\,.. . , A-i e GL„(KF k -i) 

that satisfy (13), (14), and (15) (k - 1 is substituted for m). We claim that that there exists a 
2-field F k generated over Zf°' A by at most (k - l)n 2 elements and Z E Gh„(KF k ) such that 

(Di, . . .,At-i,A0 := (Di, . . .,D k - U ZC k ) 

satisfies (13), (14), and (15) (k is substituted for r). We need to construct a 2-field F k and 
Z e GL n (KF k ) such that 

o- k (Ai)ZC k = <Pi(ZC k )A h UK?, (16) 

cr k (BdZC k = <9,(ZQ ) + ZC k B/, 1 sC / < m, (17) 

cr*( D,)ZC, = o-i(ZC k )Di, l^i^k-1. (18) 

Expanding the right-hand sides of (16), (17), and (18) using (9) and (10), we have 

<Pi(ZC k )Ai = 0i(Z)0i(Q)Ai = frWo-kiAdCk, l^i^q, 
di(ZC k ) + ZC k Bi = di(Z)C k + Z(di(C k ) + C k Bi) = di(Z)C k + Zo- k {B-)C k , 1 < i < m, 
o-i(ZC k )Di = o-i(Z)o-i(C k )D h 1 < i < * - 1. 

Therefore, (16), (17), and (18) turn into 

<t>i(Z) = o-tiAdZo-tiAi)- 1 , l^i^q, (19) 

di(Z) = [o- k (Bi),Z], Ui< m, (20) 

(Ti(Z) = (otfAJZortA) -1 ) (^(AOCtDrV^Ci)" 1 ) , 1 < i < Jt - 1. (21) 

Let us show that, for all i, 1 ^ i ^ k — 1, 

o-r 1 (^(A-)QZ)rVi(Q)- 1 ) (22) 
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satisfies (19) as an equation in Z. Indeed, for all j, 1 ^ j ^ q, we have 

<f>j (o~k(Di)C k DJ l o-i(C k y^ = (T k (4> j (D i ))(t>j(C k )(l> j (D i r X (r i ((t>j(C k )r x = 

= cr^jiDdAj^iAjr^jiC^jiDd-'cr^jiCk))- 1 = 

= ( r i (cr i (A y ))^(DO^(^)"Vi(Q)^ ! )" 1 o-i(^(Q))" 1 = 
= (r,(cr (i (A i ))^(D,)C t A;V j (A) _ V / (^(Q))" 1 = 

= cr,(^(A i ))^(D,)C,D7V ; (A y )-V ; (^(C,))" 1 = 
= a ,(cr k (A j))cr k (Di)C k D] 1 cr i(C k y l o- ,(cr k (A j))' 1 , 

which implies the claim. Moreover, we will prove that (22) satisfies (20) as an equation in Z. 
Indeed, for all j, 1 ^ j : ^ m, 

(TicT^B^akiDdCkD^aiiCt)- 1 - cr k (D i )C k DJ 1 cri(C k y 1 cri(a- k (Bj')) = 
= o- k (dj(Di) + DiB^CtDj'cniC,)- 1 - cTkiDdCkDj'cniCkr'cnicrdBj)) = 
= djicTkiDiVCkDj'cniCkr 1 + CTkiDdcTkiB^CkD^cniCkr 1 - 

- (TkiDdCkD^cn (cl l (T k {Bjj) = 

= djicTkiD^CkDj'cniCkr 1 + (TkiDdcrkiB^CkDi'tTiiCkr 1 - 

- cTkiDdCkDj'cTi {C k y dj{C k )C k l + Bf- k y ) = 

= djio-kiDiVCkD^o-iiCkT 1 + (TkiDdcrkiB^CkD^cniCkr 1 - 

-crkiDdCkDj'crjiBj^iiCkr 1 +(Tk(D i )CkDJ l dj( ( T i (Ckr 1 ) = 
= dj^kiDi^CkDj'aiiCk)- 1 + o- k (Di)(dj(C k ) + C k Bj)DJ l (Ti(C k y l - 

-o-kiDdCkD^o-iiB^cniCk)- 1 +cr k (D i )C k D; l dj(cr i (Ckr 1 ) = 
= djiniD^CkDj'criiCk)- 1 + <r k (D i )d j (C k )DJ 1 o- i (C k r 1 + 

+ (TkiDdCk (BjDJ 1 - Dfo-iiBjj) cniCkT 1 + crk(Di)CkDJ l dj (cr,(Q) _I ) = 
= dj(crk(Di))CkDJ l o-i(CkT l + crk(Di)dj(Ck)DJ l cri(Cky l + 

+ o- k (Di)C k dj (Di^cniCk)- 1 + (TkiDdCkDj'dj (^(Ck)-') = 

= dj((Tk(D i )CkD; 1 cr i (Ckr l ). 

Let us show that, for any {O U 2, A}-ring R and Z e GL„(R) satisfying ( 19) and (20), and, for 
all i, Ui<)t, 

ZcrJ 1 (crkiDdZcrkiDir'y 1 
also satisfies (19) and (20). Indeed, for all j, 1 ^ j ^ q 

a-ityjfecrj 1 (crkiDdZ-'cTkiDd- 1 ))) = 

= cr i (cr jt (A J )) C r i (Z)cr i (cr i (A J ))- 1 (r it (0 J (Z) i ))cri(A J )Z- 1 cri (ajV/A) -1 ) = 
= ( r,(^(A i ))cr,(Z)cr,((r,(A,))" 1 cr t (cr,(A i )A)Z- 1 ^ (DJ 1 cr^A j)' 1 ) = 
= cri (crk(Aj)ZcTT l (cTkiDdZ-'crkiDir^cTkiAjr 1 ) . 
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Moreover, for all /', 1 ^ j ^ m, 

[cr /(o" k(B ))), o-ffiakiDdZ^a-kiDi)- 1 ] = 

= ^(^(B^cr^D^Z-V^A)" 1 -cr ! (Z) ( r t (D,)Z-V,(D,)- 1 ( r,( f rA(B./)) = 

= (Ti(dj(Z) + ZcrdBj^iDdZ-'crdDir 1 - criiZ^DdZ' 1 ^ (DJ l dj(Di)D; 1 + B jDJ 1 ) = 

= 5 j (^(Z))^(D,)Z- 1 ( r,(A)- 1 + (TiiZp-Ao-i&dDi^otfPiT 1 - 

- cr i (Z)(Ti(Z) i )Z- 1 cr t (B J ) C r i (D i )" 1 + °" ;(Z)cr k (D i )Z~ 1 dj [cr ^(D/) -1 ) = 
= SXcrXZ))^^^,^- 1 ^^^)- 1 + a-i(Z)cTk{dj(Di) + DjA^artA) -1 - 

- crK^cr^DOZ-V^B.Ocr^A)" 1 + o" /(Z)cr k {Di)Z~ x dj (cr fc(Z);) _1 ) = 

= 5Xcr,(Z))cr,(D,)Z- 1 cr i .(A)- 1 + (Ti(Z)d j^Di))?' 1 ^{0^ + 

+ CT,(Z)cr,(D,) [cr^B,)^- 1 ] 0" fc(-Dj) _1 + CTiCZ^CAOZ- 1 ^- ((^(A)" 1 ) = 

= d j ((T i {Z)(r k {D i )Z- 1 (T k {D i r l ). 

Let us now demonstrate how find an invertible solution of (19), (20), and (21). Let V be the 
F&-vector space of matrix solutions of (19) and (20) in M^A^Fn) and let dim^, V =: d ^ n 2 . 
By the above, d ^ 1. Choose an -basis 



ZnU 



Zln,l 



Znl.d 



of V so that 



Zii.i 



Zln.l 



\Zwi,l/ \Znn,dJ 



= ^(DOC^rViCC*) -1 6 GL„(^Fn) 



(23) 



VZnl.l . ■ . Znn,lJ 

and, in this basis, (21) is of the form 



We now let 



St = {cr t , . . . , cr,-}, Ft := Ft-if*,, 1 sC s ^ d} T , cr,- 



1 < i ^ k- 1, 



where jt s , 1 ^ s ^ n, are Z^-indeterminates. We need to prove that the action of cr, , 1 ^ i ^ k— 1, 
is well-defined. For this, it is sufficient to show that, for all m and v, l^w, 



cr u (E v )E u = cr v (E u )E v . 



For this, it is sufficient to show that (21) satisfies the integrability conditions. For all u and v, 
1 ^ u, v ^ k — 1, we have 

cr u (o- k (D v )Z)cr u [c k D- l cr v {C k y l ) = 

= a k (<t v {D u )D v D- x ) cr k (D u )ZC k D- l cr u (C k r'cr u [c k D- 1 cr v (C k y l ) = 
= a x {(r k {D^)cr k (D x )ZC k D- l (T u {D x r'a- r {(T u (C k )r l = 
= cr^D^Z^C^D.D-'cr^D^crvicrJC^)- 1 = 
= a v {(T k (D u )Z) o~ v [c k D- u l (T u (C k r l ) . 

It remains to note that, since x\, . . . , Xd are transcendental over KF k -\, the matrix 



(zu,\x\ + ... + Zu,dXd 



Z := 



VZral,l*l + . . . + Z n l,dX d 

by (23) and satisfies (19), (20), and (21). 



Zln,\X\ + ...+ Z\n,d*d 



-■■mi. 

1*1 + . . . + Znn,d x d) 



e GL n (KF k ) 



□ 
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